﻿ Seria Fourier Analiza spectrală a semnalelor periodice Daca descompunem semnalul de intrare periodic intr-o serie de componente mai simple, putem calcula raspunsul la fiecare componenta si face sinteza raspunsurilor partiale In domeniul frecventa: seria Fourier http://shannon etc upt ro/teaching/ssist/Cap4 pdf 1 Răspunsul sistemelor continue liniare şi invariante în timp la exponenţiala complexă de modul unitar h(t)  0jt    0jt x t e   y t h e d           jtj          dehety00     H() Transformata Fourier a raspunsului la impuls h, 0 calculata in : depinde de  si h 00 2 1 0jt 0jt x t e        0y t e H h(t) Functie proprie a SLIT Valoare proprie a SLIT  jj      H h e d H e                  00jtjt 0        00y t e H H e h(t) jt  tj kk          kky t a H e keatx kk Daca semnalul de intrare este o combinatie liniara de exponentiale complexe  iesirea : o combinatie liniara de exponentiale complexe kkj t j t     k k ky t a S e a H e 3 kk Transformari ortogonale •Produsul scalar al vectorilor TT xy     1 2 1 2 nnx x x y y y ; * 1y * * * *2y xy       1 2 1 1 2 2, n n nx x x x y x y x y   * ny  •Produsul scalar al functiilor din L2 [a,b] b * ,x t y t x t y t dt  a 4 2 •Se observa ca se indeplinesc urmatoarele conditii: * i) , , ,x y y x ii ) x , y  z  x , y  x , z , iii) x y x y , , ,  , ,*iv) ,x y x y C      n m n m * v) , , x y x y        k k l l k l k l     1 1 1 1k l k l •Norma ||x|| este finita (spatiul L2): n 2 2 2 2 2 x x x       1 1 1, kx x x x 1k 22b x x t dt   a Un spatiu vectorial cu norma definita prin produsul scalar este un spatiu Hilbert (teoria aproximarii) 5 Vectori ortogonali (perpendiculari) •Pentru doi vectori bidimensionali x i j y i j    x x y y1 2 1 2 ; x y x y  •Produsul scalar : , cos ,xy cos   xy -unghiul dintre vectori Conditia de ortogonalitate : Produsul scalar sa fie zero ,0  x y x y6 3 Functii ortogonale •Vom considera doua semnale definite pe 2 (0,T), cu T=2/– spatiul L 000      00cos ; sinx t t y t t •Produsul scalar este: 00TT 1            0 0 0 0 0cos ,sin cos sin sin 2t t t t dt t dt  2 00 0T     0cos 21 cos 4t    0 0044 0 7 Spatiul Hilbert •Un sistem U={u} de vectori ortogonali doi cate doi se k spune ca este complet in spatiul Hilbert, H, daca nu exista nici un vector xH-U, care sa fie ortogonal pe toti vectorii din U (doar vectorul 0): , 0 0, if u x x x H U    k  •Un sistem complet U formeaza baza ortogonala in spatiul Hilbert Pentru orice element x din H, exista o dezvoltare unica de forma , x H x a u   kk k8 4 Exemple •Daca multimea elementelor din U, n, este finita: spatiu Hilbert finit dimensional, cu dimensiunea n •Daca multimea este numarabila dar infinita: spatiu Hilbert infinit dimensional •Versorii {i, j, k} formeaza o baza in spatiul tridimensional, cu n=3 jkt •Multimea functiilor {e 0}| cu frecventa k – kZ0 o baza infinit dimensionala pentru semnale periodice in timp continuu, de perioada T 0 9 Teorema lui Pitagora in spatiul Hilbert Relatia dintre distanta si produsul scalar •Fie diferenta intre doi vectori din spatiul Hilbert d x y  •“Distanta” dintre ei: 2 2 x y x y ,d   •Avem in general: 2 2 ,d x y x y   22  2 Re ,x x y y    22 2 ,d x y x y  •Daca x si y sunt ortogonali  10 5 2 Exemple, L •Norma pentru semnalele ortogonale     00cos si sintt 00TT  21 cos 211tTt  TT 00020             0 0 0 0cos sin 2x t t dt dt t  0002 2 2 2 2 •Distanta dintre semnale este (cf teoremei lui Pitagora) 2TT 00      0 0 0cos ,sin d t t T 22 •Semnalele ce nu sunt ortogonale nu satisfac teorema lui     Pitagora Ex:  00cos si costt 0T 2        0 0 0 0cos , cos cos / 2 t t tdt T  0 2 11       0 0 0 0 0 0cos , cos / 2 / 2 2 / 2 2 d t t T T T T cos ,sin cos , cos d t t d t t     0 0 0 0 12 6 Inegalitatea lui Schwartz in spatiul Hilbert (Cauchy- Bunyakovsky-Schwarz) Egalitatea are loc daca si numai daca x si y sunt ,x y x y dependente liniar, y=kx, unde k=const •Semnale ortogonale L2                 00cos , sin , 0x t t y t t x t y t •Produsul normelor este: T T T 0 0 0 x t y t        2 2 2 •Inegalitate 0 0 si ω<0, (k, ) 0k  sgnk    k 2 Functie impara 34 Seria Fourier exponentiala 17 Alte forme ale relatiei lui Parseval •Seria Fourier exponentiala :  22 2 21    2P x t dt c c c   0kk  T     000kkT •Forma trigonometrica si armonica 2 2 2 2a b A 221k k k          00P a x t dt A T 0112 2 2kkT 0 •Exemplu Puterea semnalului rectangular: 0/2T 211 42P x t dt dt     4T35 000T Diagrama spectrala de putere folosind 2 seria Fourier armonica (k, A/2) 0k Semnalul rectangular Frecvente exclusiv pozitive Se recomanda folosirea unei reprezentari logaritmice pentru putere ce avantajeaza reprezentarea puterilor mai putin semnificative 36 18 Diagrama spectrala de putere folosind 2 seria Fourier exponentiala (k, |c|) 0k Frecvente pozitive si negative 37 •Pentru semnale de banda nelimitata : –Banda de frecventa este infinita –Puterea scade cu cresterea frecventei, tinde spre zero pt frecvente ce tind la infinit •Banda efectiva de frecvente = gama pozitiva de frecvente ce contin un procentaj semnificativ al puterii semnalului •In acest caz, in banda 9 se gaseste 96,5% din 0 puterea semnalului 38 19 Fenomenul Gibbs •Fizicianul Albert Michelson a construit un analizor de spectru in 1898 •La iesirea filtrelor analizorului, a masurat amplitudinile componentelor spectrale, conform teoriei •Cand a incercat sa recompuna prin insumare, semnalul initial, a observat ca apare o “problema”, la semnal anume L-a rugat pe Gibbs sa ii explice acest fenomen 39 Semnalul analizat (de banda nelimitata): rectangular cu factor de umplere 0 5, fara componenta continua 4 1 1         0 0 0sin sin 3 sin 5 x t t t t  35  Printr-o trunchiere in frecventa, pastrand primele n armonici, de ordin impar, semnalul este aproximat cu unul de banda limitata: 4 1 1 1           0 0 0 0sin sin 3 sin 5 sin 2 1x t t t t n t  3 5 2 1n  40 20 02nt 2 sin 2u       0Si 2x t du n t  0u •Si(x) – sinus integral, functie impara http://mathworld wolfram com/SineIntegral html /2 x sinu Si ; Si Six du x x           0u -/2  limSix  41 2x Fenomenul Gibbs •Gibbs a aratat ca trunchiind semnalul rectangular cu factor de umplere 0 5, si pastrand n armonici de ordin impar, 4 1 1 1           0 0 0 0sin sin 3 sin 5 sin 2 1x t t t t n t   3 5 2 1n  •Se obtine 2nt 0 2 sin 2u    0Si 2x t du n t  0u •Semnalele de banda nelimitata nu pot fi perfect aproximate cu semnale de banda limitata 42 21 O unda rectangulara cu T=1s, cu 2nf=80f 000 Asimptote orizontale: 1, -1 43 •Prima supracrestere (maximul oscilatiei), de 1 18 V apare la momentul t=6,25ms << T m0 •Oscilatii in vecinatatea zonei de tranzitie •Depasirea valorii 1 este de 9% •Timpul de ridicare (rise time) 21 2tt   rm MMf ω=(2n-1)ωaproximativ 2nω– cea mai mare M0 0 frecventa retinuta 44 22 Semnale trunchiate pentru 21 si respectiv 45 armonici Se poate observa ca oscilatiile nu scad ca si amplitudine, dar frecventa lor creste Semnalul de aproximare converge 45 in medie patratica catre semnalul x(t) Distributia Dirac periodica  1            Tkt t kT c   0kT 2 22TT 1 1 1jk t T       kTc t e dt t dt T T T 22TT •Pentru [-T/2,T/2] , (t)= (t) T  01jk t      Tt t kT e T   kk46 23 Proprietatile seriei Fourier exponentiale •Coeficientii seriei Fourier a semnalului x, de perioada T x       kx t c 01jk t   kc x t e dt  TT •Descompunerea Fourier  j k  t 0 a p t (aproape peste tot)x t c e   k k   47 1 Liniaritatea •semnalele x(t) si y(t) periodice cu perioada T : xy  x t c y t c ,            kk xy           kkax t by t ac bc 2 Deplasarea în timp 00jk tx        0kx t t e c          000 0 011jk tjk t jk tx           0kkc x t t e dt x e d e c  TTTT •Deplasarea in timp  modulatie cu exponentiala complexa 48 24 3 Conjugarea complexă •Conjugarea complexa in timp  reflectarea in domeniul frecventa si conjugarea complexa * *x      kx t c * *11      00*j k tjk t            kkc x t e dt x t e dt c  TT 0000TT  4 Reflectarea semnalului     0011jkjk tx            kkc x t e dt x e d c  TTTT x            kx t x t c Reflectarea in timp  reflectare in frecventa •49 5 Scalarea variabilei timp •x(t) - perioada T  x(at), perioada T/|a|    012jk t        00;kc x at e dt a  //TaT T a  01xjk       kkc x e d c  TT x       kx at c 6 Modularea semnalului •Modulatia in timp  deplasare in domeniul frecventa      000 0 011kkjtjk t jk tx c x t e e dt x t e dt c     0kkk TTTT jk tx50 00 x t e c   0kk 25 Dualitatea timp-frecventa •O operatie in timp  alta operatie in frecventa: –De exemplu: modulatie in timp  deplasare in frecventa •A doua operatie in timp  prima operatie in frecventa –Deplasare in timp  modulatie in frecventa •Acest comportament este numit dualitate •Reflectarea este o operatie auto-duala 51 7 Produsul a două semnale •Convolutia coeficientilor  xyxy              k n nkkx t y t c c c c n    8 Convoluţia periodică a semnalelor •Semnalele periodice nu au energie finita, si convolutia nu se poate defini Se foloseste convolutia circulara sau periodica, definita pe o perioada xy         kkz t x y t d x t y t Tc c  T Operatii duale: inmultirea ↔ convolutia •52 26 Convoluţia periodică a doua semnale rectangulare, cu factor de umplere diferit •Efect de circularitate 53 9 Derivarea semnalului •Dupa diferentiere, componenta continua=0 Semnalul ramane periodic Derivarea in timp  inmultirea spectrului cu jkω 0   xdx t jk c     0k dt 10 Integrarea semnalului •Pentru ca semnalul sa ramana periodic dupa integrare, componenta continua trebuie sa fie nula Integrarea in timp  inmultirea spectrului cu 1/jkω 0 tx xkc     0 0x d c  jk  054 27 11 Semnale reale Seriile componentelor para si impara *  •x(t) semnal real; kkcc •Componentele para x(t) si impara x(t) pi •Spectrul componentei pare x(t) –real p •Spectrul componentei impare (semnal real) x(t) – i pur imaginar     xx t x t   Rex t c     kp  2  xx t x t Imx t j c    ki55 2 28 